Harmonic sensor

ABSTRACT

A system including a weighting function module, a frequency sensor module, and a control module. The weighting function module is configured to receive a plurality of samples of data, the data having been obtained based on a head having read the data from a storage medium while the head was at a given height above the storage medium, and to apply a weighting function to the plurality of samples to generate weighted samples. The frequency sensor module is configured to estimate, based on the weighted samples, (i) a first magnitude of a first frequency component, and (ii) a second magnitude of a second frequency component. The control module is configured to estimate, based on the first magnitude and the second magnitude, the height of the head above the storage medium.

CROSS-REFERENCE TO RELATED APPLICATIONS

This is a continuation of U.S. patent application Ser. No. 13/442,033(now U.S. Pat. No. 8,422,160), filed on Apr. 9, 2012, which is acontinuation of U.S. patent application Ser. No. 13/190,036 (now U.S.Pat. No. 8,154,820), filed on Jul. 25, 2011, which is a continuation ofSer. No. 12/581,262 (now U.S. Pat. No. 7,986,487), filed on Oct. 19,2009, which claims the benefit of U.S. Provisional Application No.61/110,331, filed on Oct. 31, 2008. The entire disclosures of the aboveapplications are incorporated herein by reference.

FIELD

The present disclosure relates to a data storage medium. Moreparticularly, the present disclosure relates to using a harmonic sensorto determine fly height of a read/write head over the data storagemedium.

BACKGROUND

The background description provided herein is for the purpose ofgenerally presenting the context of the disclosure. Work of thepresently named inventors, to the extent the work is described in thisbackground section, as well as aspects of the description that may nototherwise qualify as prior art at the time of filing, are neitherexpressly nor impliedly admitted as prior art against the presentdisclosure.

Complex signals can be represented in a time and/or a frequency domain.Signals may be transformed from the time domain to the frequency domainand vice versa. Transformation from the time domain to the frequencydomain may be implemented by a Fourier transform. The Fourier transformprovides a spectrum of frequency components (i.e., spectra), referred toas “harmonics,” each of which has a corresponding magnitude. Consideredtogether, the signal spectra constitute a frequency domainrepresentation of a given signal.

One notable use of such a transformation, is to electronically determinethe fly height of a read/write head in a hard disk drive (HDD). The flyheight is the height of the head relative to a recording medium.

Methods for detecting fly height generally employ detecting a knownpattern written on the storage medium. In response to detecting thepattern, the head produces a signal. A ratio in harmonic magnitude ofthe fundamental frequency (i.e. first harmonic (Fc)) and correspondingthird harmonic (3Fc) of the signal may directly correlate to a change infly height (ΔD) of the head. This relationship may be expressed as:ΔD=|3Fc|÷|Fc|

In storage media systems, fly height of the head is often maintainedwithin a prescribed region. If the fly height is too great, the signalstransmitted and read by the head may be too weak for accurate datastorage and retrieval. If the fly height is too small, the head mayphysically contact the medium, causing damage to the medium and/or lossof the data stored on the medium.

SUMMARY

A fly height control system includes a weighting function module. Theweighting function module applies coefficients from a window function ina time domain to digital samples to provide weighted samples. A first ofthe coefficients is applied to a first of the digital samples and asecond of the coefficients is applied to a second of the digitalsamples. The first of the coefficients has a first magnitude that isdifferent than a second magnitude of the second of the coefficients. Thesystem also includes a harmonic sensor module to estimate a magnitude ofa frequency in a frequency domain based on the weighted samples. Thesystem also includes a head height control module determine a height ofa read/write head over a storage medium based on the magnitude of thefrequency.

In other features, the frequency in the frequency domain is a firstfrequency; and the harmonic sensor module is a first harmonic sensormodule. The system further comprises a second harmonic sensor module toestimate a magnitude of a second frequency in the frequency domain basedon the weighted samples. The head height control module determines theheight of the read/write head over the storage medium based on themagnitude of the first frequency and the magnitude of the secondfrequency.

In other features, the first harmonic sensor module performs a firstdiscrete Fourier transform on the weighted samples to estimate themagnitude of the first frequency. The second harmonic sensor moduleperforms a second discrete Fourier transform on the weighted samples toestimate the magnitude of the second frequency.

In other features, the head height control module determines the heightof the read/write head over the storage medium based on a ratio of themagnitude of the first frequency to the magnitude of the secondfrequency. The head height control module adjusts the height of theread/write head over the storage medium when the ratio indicates thatthe height of the read/write head over the storage medium is outside ofa predetermined range. The first frequency corresponds to a firstharmonic of at least one data pulse from the read/write head, and thesecond frequency corresponds to a third harmonic of the at least onedata pulse.

In still other features, the systems and methods described above areimplemented by a computer program executed by one or more processors.The computer program can reside on a computer readable medium such asbut not limited to memory, nonvolatile data storage, and/or othersuitable tangible storage mediums.

Further areas of applicability of the present disclosure will becomeapparent from the detailed description, the claims and the drawings. Thedetailed description and specific examples are intended for purposes ofillustration only and are not intended to limit the scope of thedisclosure.

BRIEF DESCRIPTION OF DRAWINGS

The present disclosure will become more fully understood from thedetailed description and the accompanying drawings, wherein:

FIG. 1 is a high level schematic diagram of a read/write head for a harddisk drive;

FIG. 2 illustrates a plot of magnitude versus frequency for variouspulse densities;

FIG. 3A is a block diagram of a hard disk drive system according to thepresent disclosure;

FIG. 3B is a block diagram of a R/W channel for the hard disk drivesystem according to the present disclosure;

FIG. 3C is a block diagram of a R/W channel for the hard disk drivesystem according to the present disclosure;

FIG. 4 illustrates examples of time domain window functions according tothe present disclosure;

FIG. 5 illustrates a comparison of window function spectra according tothe present disclosure;

FIG. 6 illustrates harmonic ratios for a 4T pattern according to theprior art;

FIG. 7 illustrates magnitude versus frequency for random patterns absentwindow function weighting;

FIG. 8 illustrates a magnitude of the output of a first harmonic sensormodule divided by a magnitude of the output of a second harmonic sensormodule for various patterns absent window function weighting; and

FIG. 9 is a block diagram of a head height control module according tothe present disclosure.

DESCRIPTION

The following description is merely exemplary in nature and is in no wayintended to limit the disclosure, its application, or uses. For purposesof clarity, the same reference numbers will be used in the drawings toidentify similar elements. As used herein, the phrase at least one of A,B, and C should be construed to mean a logical (A or B or C), using anon-exclusive logical OR. It should be understood that steps within amethod may be executed in different order without altering theprinciples of the present disclosure.

As used herein, the term module may refer to, be part of, or include anApplication Specific Integrated Circuit (ASIC), an electronic circuit, aprocessor (shared, dedicated, or group) and/or memory (shared,dedicated, or group) that execute one or more software or firmwareprograms, a combinational logic circuit, and/or other suitablecomponents that provide the described functionality.

Typically, predetermined test patterns are written and then read-out ofa storage medium as pulses. The pulses are used determine and adjust flyheight of a transducer over the storage medium. An example of atransducer is a read/write head of a hard disk drive (HDD).

A read/write (R/W) channel (also referred to as a recording channel) ofthe present disclosure includes harmonic sensors. The harmonic sensorsreceive analog to digital samples from an analog to digital converter(ADC) that are based on pulses from the head. The harmonic sensorsprovide outputs that are used to estimate fly height based on random orpseudo-random data.

In one embodiment, the R/W channel applies a time-domain window over theADC samples before they are received by the harmonic sensors. Thetime-domain window weights the ADC samples according to a predeterminedpattern. Weighting may be provided by multiplying coefficients providedby the window function to incoming samples. The harmonic sensors use thewindowed samples to estimate spectrums of signals from the head. Thewindow reduces dependency of spectrum estimation of the harmonic sensorson the particular random data pattern written. An example of a randompattern is a data pattern written to a data field.

The harmonic sensors may provide discrete Fourier transforms (DFTs) atvarious frequencies based on the time-domain window and ADC samples. Ahead height control module may sum the DFTs for the various frequenciesand estimate and adjust head height based on the sums, as will bediscussed herein.

Referring now to FIG. 1, when writing data, a write current is generatedthat flows through the head 106. The write current is switched toproduce a magnetic field having a positive or negative polarity. Thepositive or negative polarity is stored by the medium 107 as a magnetictransition representing data. In other words, a transition 130 indicatesa transition from a positive polarity to a negative polarity or viceversa. When reading data, the head 106 provides a pulse based on thetransition 130.

The pulse may correspond to a single transition pulse or a dibittransition pulse. For a dibit pulse, a binary “1” may be encoded as apair of magnetic transitions spaced one bit apart. A binary “0” may beencoded as the absence of any transitions in a particular location onthe medium 107.

In FIG. 1, the head 106 is illustrated at two different heights relativeto the medium 107. The first height D1 corresponds to the head 106flying close to the medium 107. The second height D2 corresponds to thehead 106 flying further away from the medium 107 than at the firstheight.

As the head 106 flies close to the medium 107, magnetic field lines 134are tightly bunched causing a sharp transition 130. In other words, thehead 106 takes a relatively short amount of time to pass the field lines134 over the transition 130 and provide a pulse based on the transition130. Therefore, the pulse width (PW) may be relatively small.

When the fly height is higher, the field lines 136 are further apart,and the transition 130 maybe more gradual. In other words, the head 106takes a relatively large amount of time to pass the field lines 136 overthe transition 130 and provide a pulse based on the transition 130. Thepulse width for the gradual transition may be larger than the pulsewidth for the sharp transition. Increases in pulse width havecorresponding increases in pulse density. Thus, the pulse density, alsoreferred to as channel bit density (CBD), may be higher for the gradualtransition. Fly height may therefore be determined based on pulsedensity.

In one embodiment,CBD=PW₅₀ /T  (1)T is the symbol duration and PW₅₀ is the pulse width at thehalf-amplitude point of the pulse. PW₅₀ for the head 106 may beapproximated by:PW₅₀≈√{square root over (2g ²+4(d+a)²)}  (2)where g is the gap of the head, a is the transition length, and d is thefly height. Once pulse density is known, equations (1) and (2) can beused to determine the fly height d. Thus, the higher fly heightcorresponds to a higher pulse density.

When the head 106 encounters a transition on the medium 107, the head106 generates a signal s(t) that can be modeled by:

$\begin{matrix}{{{s(t)} = \frac{1}{1 + \left( \frac{2\; t}{P\; W_{50}} \right)^{2}}},} & (3)\end{matrix}$where t is a time variable. The signal s(t) can be written as:

$\begin{matrix}{{y(t)} = {{\sum\limits_{m}{\left( {u_{m} - u_{m - 1}} \right){s\left( {t - {mT}} \right)}}} + {{n(t)}.}}} & (4)\end{matrix}$U_(m) is the data stream written onto the medium 107, spaced by time T,and n(t) is the noise in the signal. After ADC sampling, thediscrete-time signal (y_(k)) can be written as:

$\begin{matrix}{y_{k} = {{y\left( {k\; T} \right)} = {{\sum\limits_{m}{\left( {u_{m} - u_{m - 1}} \right){s\left( {{k\; T} - {m\; T}} \right)}}} + {{n\left( {k\; T} \right)}.}}}} & (5)\end{matrix}$Y_(k) depends on s(t), which is a function of fly height d.

Pulses from the head 106 may be translated by harmonic sensors from thetime domain to the frequency domain. A ratio of a low frequency in thepulse to a high frequency in the pulse, as represented in the frequencydomain, may indicate a corresponding magnitude of pulse density. Forexample, a low frequency harmonic (e.g., first harmonic) may beselected, and a higher frequency harmonic (e.g., third harmonic) may beselected. Therefore, the ratio of harmonics may also be used todetermine fly height, as will be further discussed regarding FIG. 2.

Referring now to FIG. 2, as pulse density decreases, higher frequencieshave more energy, and the ratio of the first harmonic to the thirdharmonic goes down. FIG. 2 illustrates a plot of magnitude versusfrequency for various densities. Magnitudes can be calculated from a sumof Fourier series coefficients output from the harmonic sensors, as willbe discussed herein.

Table 1 (provided below), which corresponds to FIG. 2, indicates that asthe pulse density (K) increases, the ratio of the first harmonic (e.g.,⅛ Fs) to the third harmonic (e.g., ⅜ Fs) also increases.

TABLE 1 K f = 1/8 Fs f = 3/8 Fs Ratio 0.80 0.87 0.30 2.96 1.00 0.82 0.174.87 1.20 0.76 0.09 8.77A first density for the pulse is represented by curve 150, a seconddensity for the pulse is represented by curve 152 and a third densityfor the pulse is represented by curve 154. The first density is higherthan the second density, which is higher than the third density. Aspulse density decreases, higher frequencies have more energy, and theratio of the first harmonic to the third harmonic goes down. In otherwords, the density may be derived from the ratio of the first harmonicto the third harmonic.

In one embodiment, the distance between the head 106 and the medium 107increases when this occurs, the amplitude of a signal, and particularly,a third harmonic of a pulse from the head 106 is decreased in the timedomain. Accordingly, the magnitude of the third harmonic output from theharmonic sensor is also decreased in the frequency domain.

Therefore, the spectrum of the pulse, and in particular, the roll-off inthe frequency domain, may indicate the density of the pulse. Roll-offcorresponds to the rate of loss or attenuation of a signal beyond acertain frequency. The higher the density, the faster the roll-off, andthe less energy left over in the pulse signal at higher frequencies.Therefore, fly height may also be determined based on roll-off. Thewindow functions of the present disclosure may minimize effects ofundesirable frequencies in the pulses. The window functions maytherefore allow the harmonic sensors to more accurately translate theharmonics into the frequency domain.

The teachings of the disclosure can be implemented in a hard disk drive(HDD) 400, an example of which is illustrated in FIG. 3A. The HDD 400includes a hard disk assembly (HDA) 401 and a HDD PCB 402. The HDA 401may include a storage medium 107, such as one or more platters thatstore data, and a read/write head 106. The head 106 may be arranged onan actuator arm 405 and may read and write data on the medium 107. TheHDA 401 includes a spindle motor 406 that rotates the medium 107 and avoice-coil motor (VCM) 407 that actuates the actuator arm 405. The HDA401 may also include a head height adjustment motor 410 that moves thehead 106 closer to or further from the medium 107.

A preamplifier device 412 amplifies signals generated by the head 106during read operations, and provides signals to the read/write head 106during write operations. The head 106 transfers data to and/or fromstorage medium 107 by writing transition pulses to or reading transitionpulses from the medium 107.

The HDD PCB 402 includes a R/W channel module 414, a head height controlmodule 415, a hard disk controller (HDC) module 416, a buffer 418,nonvolatile memory 420, a processor 422, a spindle/VCM driver module 424and a height driver module 426. The R/W channel module 414 processesdata received from and transmitted to the preamplifier device 412. TheHDC module 416 controls components of the HDA 401 and communicates withan external device (not shown) via an I/O interface 430. The externaldevice may include a computer, a multimedia device, a mobile computingdevice, etc. (not shown). The I/O interface 430 may include wirelineand/or wireless communication links.

The HDC module 416 may receive data from the HDA 401, the R/W channelmodule 414, the head height control module 415, the buffer 418,nonvolatile memory 420, the processor 422, the spindle/VCM driver module424, and/or the I/O interface 430. The processor 422 may process thedata, including encoding, decoding, filtering, and/or formatting. Theprocessed data may be output to the HDA 401, the R/W channel module 414,the buffer 418, nonvolatile memory 420, the processor 422, thespindle/VCM driver module 424, and/or the I/O interface 430.

The HDC module 416 may use the buffer 418 and/or nonvolatile memory 420to store data related to the control and operation of the HDD 400. Thebuffer 418 may include DRAM, SDRAM, etc. The nonvolatile memory 420 mayinclude flash memory (including NAND and NOR flash memory), phase changememory, magnetic RAM, or multi-state memory, in which each memory cellhas more than two states. The spindle/VCM driver module 424 controls thespindle motor 406 and the VCM 407. The height driver module 426 controlsthe head height adjustment motor 410 to move the head 106 closer to orfurther from the medium 107. The HDD PCB 402 includes a power supply 417that provides power to the components of the HDD 400.

Referring now to FIG. 3B, the R/W channel module 414 is shown inaccordance with one embodiment of the present disclosure. The R/Wchannel module 414 includes an analog front end (AFE) module 520, ananalog to digital converter (ADC) module 522, a finite impulse response(FIR) filter 524, a Viterbi decoder (VTB) 526, a weighting functionmodule 527 and harmonic sensor circuit (HSC) 530-1, 530-2.

In operation, the R/W channel module 414 processes the signals from thehead 106. To this end, the analog front end module 520 receives the datato/from the head 106 in the analog domain. For example, the analog frontend module 520 receives and processes pulses from the pre-amplifierdevice 412. The analog front end module 520 may include a high passfilter (HPF) 540, a variable gain amplifier (VGA) 542 and a continuoustime filter (CTF) 544.

Signals from the head 106 may be filtered and amplified through the HPF540, VGA 542 and CTF 544. In one embodiment, the VGA 542 provides gainto signals provided from the head 106. The variable gain is controlledby an automatic gain control (AGC) loop (not shown). The VGA 542provides control of the signal level such that the analog to digitalconversion of the R/W channel module 414 is optimized. For example,increases in gain may avoid quantization noise which may impact biterror rate (BER) performance.

The CTF 544 may minimize aliasing, which may occur during analog todigital conversion. In one embodiment, the CTF 544 is a multiple polelow pass filter. The filter may be target matched with the FIR filter524.

The ADC module 522, FIR filter 524 and Viterbi decoder 526 may performdata detection processes such as sequence detecting, bit recovery,parity processing, descrambling and the like. While FIG. 3B shows a FIRfilter 524 and a Viterbi decoder 526, different types of filtersdecoders and/or various other digital backend components may be used.

The ADC module 522 converts the analog signals of the R/W channel module414 to digital samples quantized in time and amplitude. The ADC module522 may be clocked by a phase locked loop (not shown), which tracks thechannel rate clock frequency.

In one embodiment, the FIR filter 524 includes a plurality of taps (orregisters) and equalizes the samples from the ADC module 522. TheViterbi decoder 526 decodes the signal output from the FIR filter 524.

The HSC 530-1, 530-2 may each include a respective bank of correlators,which implement discrete Fourier transforms (DFTs) at selectedfrequencies.

In one embodiment, the HSC 530-1, 530-2 may include DFTs that translatethe digital output of the ADC module 522 to sums of signal spectra. Forexample, the HSC 530-1, 530-2 may provide estimates of the sine andcosine components of the ADC module outputs, FIR filter outputs and/orViterbi decoder outputs according to:

$\begin{matrix}{Y_{{co}s} = {\sum\limits_{n}{{D(n)}{\cos\left( {2\pi\;{nf}_{HSC}} \right)}}}} & (6)\end{matrix}$

$\begin{matrix}{Y_{\sin} = {\sum\limits_{n}{{D(n)}{{\sin\left( {2\pi\;{nf}_{HSC}} \right)}.}}}} & (7)\end{matrix}$D(n) corresponds to ADC module outputs, FIR filter outputs and/orViterbi decoder outputs (selectable via multiplexer 670), n is thesample index, f_(FSC) is a selectable DFT frequency. In one embodiment,f_(HSC) may be selected between 1/32 Fs and ½ Fs. Each HSC 530-1, 530-2may have an independently selectable frequency.

In other words, the HSC 530-1, 530-2 may correspond to coarse spectrumanalyzers where frequencies of interest (e.g., first and thirdharmonics) are selected. The frequencies are thus selected to estimatethe energy in each frequency and to obtain a coarse picture of thespectrum signal.

In one embodiment, the HSC 530-1, 530-2 are illustrated separately forthe purpose of explaining the derivation of first and third harmonicmagnitudes. Various other HSC modules may also be included to obtainvarious other harmonic magnitudes. The HSC 530-1, 530-2 may not belimited to separate stages and can be implemented by a single HSC.

In one embodiment, a multiplexer (MUX) 670 is controlled by the headheight control module 415 to provide the HSC 530-1, 530-2 with outputsfrom the ADC 522. The multiplexer 670 is also controlled to provide theHSC 530-1, 530-2 with outputs from the Viterbi decoder 526. Outputs ofthe Viterbi decoder 526 may be referred to as non-return-to-zero (NRZ)bits. Further, HSC 530-1, 530-2 may selectively receive outputs of FIRfilter 524 based on MUX selection.

The head height control module 415 may align the estimates of the firstand third harmonics before Viterbi processing with estimates made afterViterbi processing. In other words, the head height control module 415may remove time offset between HSC estimations for ADC module outputsand Viterbi decoder outputs.

Referring now to FIG. 3C, in one embodiment, the HSC 530-1, 530-2estimate the first and third harmonics from ADC samples. HSC 530-3,530-4 estimate the first and third harmonics, respectively for thesamples after they are processed through the Viterbi decoder 526. Thehead height control module 415 may align the estimates of the firstharmonic from the HSC 530-1 and the HSC 530-3. The head height controlmodule 415 may also align the estimates of the third harmonic from theHSC 530-2 and the HSC 530-4. In other words, the head height controlmodule 415 may remove time offset between ADC module outputs and Viterbidecoder outputs.

In one embodiment, the signal spectra are the first and third harmonicmagnitudes of an applied signal. The output of HSC 530-1 is shown as(cos 1) and (sin 1). The output of HSC 530-2 is shown as (cos 2) and(sin 2). The head height control module 415 may express the magnitudesof the first and third harmonics:Fc Magnitude=√{square root over (sin₁ ²+cos₁ ²)}  (8)3Fc Magnitude=√{square root over (sin₂ ²+cos₂ ²)}  (9)

While shown as (sin) and (cos) output components, the magnitudes may becalculated directly by HSC 530-1 or HSC 530-2.

The head height control module 415 may divide the Fourier transform ofthe ADC samples by the Fourier transform of the Viterbi decoder for eachharmonic. The resulting quotient may equal an estimate of the channelresponse including an estimate of the original pulse convoluted with theoutput of the analog front end 520 (Pulse*AFE). For example,

$\begin{matrix}{\frac{{FT}({ADC})}{{FT}({NRZ})} = {{{FT}\left( {{{Pulse}\;}^{*}{AFE}} \right)}.}} & (10)\end{matrix}$

After the spectrum of the bits (i.e., FT(NRZ)) is estimated based on theoutput of the Viterbi decoder 526, the estimate is divided out from thespectrum of the samples (i.e., FT(ADC)). The head height control module415 normalizes (i.e., divides) the ADC spectrum by the NRZ spectrum togenerate a relatively smooth estimate of the channel response.

In the present disclosure, the HSC 530-1, 530-2 generate estimates basedon data that is not necessarily based on a predetermined pattern andthat may or may not be random. For example, the ADC module 522 mayprovide digitized data that corresponds to written user data written indata sectors on the medium 107. In other words, the data may correspondto normal read data and not test track data and may be provided inreal-time as the head 106 is reading the data.

The weighting function module 527 applies a weighting function in thetime domain to inputs to the HSC 530-1, 530-2. Windowing corresponds tomultiplication in the time domain, which is convolution in the frequencydomain. In the frequency domain, the spectrum is a series of narrow syncfunctions over the frequencies of interest with side lobes at nearbyfrequencies.

Referring now to FIG. 4, examples of time domain window functions areillustrated. The window functions include a triangle window 550, aflat-top window 552 and a rectangular window 554. The windows 550, 552,554 are illustrated over a plurality of samples (i.e., 4000 samples). Inpractice, the windows 550, 552, 554 may provide coefficients that aremultiplied by respective ones of the incoming samples, thereby weightingthe incoming samples.

The rectangular window 554 corresponds to a window, where all samplesare given the same weight. The flat-top window 552 weights samples thatare determined to correspond to frequencies of interest, such as firstand third harmonics, according to a sine or cosine functions. Theflat-top window 552 weights remaining samples substantially less (e.g.,0) than the samples for the frequencies of interest. The triangle window550 may weight the samples according to a linearly increasing functionand then a linearly decreasing function. The two functions may meet atsamples that are determined to correspond to a frequency of interest. Inother embodiments, almost any window function that weights samples forareas of interest greater than other samples may be used according tothe present disclosure.

Referring now to FIG. 5, a comparison of window function spectra isillustrated. The spectrum 560 of the triangle window 550, the spectrum562 of the flat-top window 552 and the spectrum 564 of the rectangularwindow 554 are illustrated. The main lobe 566 of the triangle window 550has a smaller normalized frequency range than the main lobe 568 of theflat-top window 552 and the main lobe 570 of the rectangular window 554.The side lobes 582 of the triangle window 550 have smaller magnitudesthan the side lobes 584 of the flat-top window 552 but larger magnitudesthan the side lobes 586 of the rectangular window 554. In oneembodiment, the frequency range of the main lobe and the magnitude ofthe side lobes may be minimized to minimize influence from frequenciesthat are not of interest.

In one embodiment, as in FIGS. 3B-C, the weighting function module 527includes an integrator module 790 to implement the window function 550.For example, the integrator module 790 may provide a linearly changingoutput based on a constant. The integrator module 790 may thereforeintegrate the constant until the mid-point (e.g., halfway through thetotal amount of samples for a pulse) of the window 550 is reached. Asubtraction module 792 may then subtract out integral values until zero(or some other initial state) is reached. The weighting function module527 may multiply the changing coefficients by incoming samples.

Referring now to FIG. 6, as mentioned, dedicated tracks were, in thepast, typically written on media with repeating test patterns. Therepeating test patterns were then read out to determine density. Anexample of a repeating pattern is a 4T pattern, which includes four onesand four zeros repeating. T is the duration of one bit.

Previous HSC modules would generate estimates based on the 4T pattern asoutput from an ADC module. An example of harmonic ratios for HSC moduleoutputs for various pulse densities is shown in Table 2 (provided below)and FIG. 6. In FIG. 6, HSC module harmonic ratios 600, 602, 604 areillustrated for a 4T pattern for densities K=0.8, 1 and 1.2,respectively.

TABLE 2 K f = 1/8 Fs f = 3/8 Fs Ratio 0.80 6.76 1.557 4.34 1.00 6.3590.915 6.95 1.20 5.902 0.471 12.53

Referring now to FIG. 7, the HSC 530-1, 530-2 may provide estimates thatare not focused on frequencies of interest absent weighting of pulsedata. For example, magnitude versus frequency for three random patterns650, 652, 654 are illustrated. As can be seen, it may be difficult torender accurate approximations of pulse density based on the randompatterns as compared with the 4T patterns of FIG. 6. In other words,FIG. 7 provides a graph of estimates for a system that does not includethe 4T patterns of FIG. 6 or the weighting function module 527 of thepresent disclosure. The addition of the weighting of pulse data by theweighting function module 527 clarifies frequencies of interest, incontrast to the random patterns of FIG. 7. This is because particularfrequencies are weighted by the weighting function module 527 so thatother non-weighted frequencies do not substantially affect frequencyanalysis, unlike in FIG. 7.

Referring now to FIG. 8, for various random patterns, the magnitude ofthe outputs 700 of HSC 530-1 may be divided by the magnitude of theoutputs 702 of HSC 530-2. The division results in ratios 704. FIG. 8illustrates one hundred random patterns where each pattern has a lengthof one typical data sector. Each circle in FIG. 8 represents ameasurement of the magnitude of the first harmonic, the third harmonic,or a ratio of the first and the third harmonics.

The ratios 704 include several outliers 708-1, 708-2, . . . , and 708-N(referred to as outliers 708). Further, from pattern to pattern, thereare variations. The outliers 708 and variations are caused by changes inthe pattern of a pulse from trial to trial. Thus, the amount of energyin the frequencies of interest (such as the first and third harmonics)changes for particular areas of a disk due to random data. The weightingfunction module 527 is included in the present disclosure to smooth outthe curve 705 corresponding to the ratios 704 and remove the outliers708.

In this manner, the reconstructed spectrum from random data may be usedto estimate harmonic ratios (i.e., ratios of the magnitude of the outputof the HSC 530-1 to the magnitude of the output of HSC 530-2, asrepresented by curve 705). As mentioned, a ratio of a low frequency inthe pulse to a high frequency in the pulse, as represented in thefrequency domain, may indicate a corresponding magnitude of pulsedensity, which corresponds to head height.

Referring now to FIG. 9, the head height control module 415 controls thedistance between the head 106 and the storage medium 107. A head heightestimation module 800 estimates current head height based on windowadjusted signals from the HSC 530-1, 530-2, such as those illustrated incurve 705 of FIG. 8. For example, the head height estimation module 800determines pulse density based on the HSC module signals and comparesthe pulse density to a previous pulse density. The head heightestimation module 800 may determine changes in height based ondifferences between the stored pulse density and the current pulsedensity based on the HSC signals.

A head height adjustor module 802 determines what, if any, head heightadjustments should be made based on the head height estimation from thehead height estimation module 800. The head height adjustor module 802may compare the estimated head height to a predetermined range ofheights. If the fly height is too great (i.e., above the predeterminedrange), the signals transmitted and read by the head may be too weak foraccurate data storage and retrieval. If the fly height is too small(i.e., below the predetermined range), the head may physically contactthe medium. Contacting the medium may cause physical damage to themedium and/or loss of the data stored on the medium. For example, if thehead height adjustor module 802 determines that the head 106 is flyingtoo low to the medium 107 such that impact may be imminent, the headheight adjustor module may adjust the height of the head 106 upwards.

The head height control module 415 may also include a selector module804 that selects frequencies for the HSC 530-1, 530-2 and that controlsthe multiplexer 670.

In one embodiment, the head height control module 415 may determine theratio of the first harmonic to the third harmonic from HSCs and maytranslate the ratio into a corresponding pulse density. In analternative embodiment, a single HSC estimates the sine and cosinecomponents at one frequency. Head height is then measured by comparingthe single frequency response to a pre-calculated or pre-calibratedlevel. Relative changes in head height may also be detected by detectingrelative changes in the response at this single frequency.

Any processes descriptions or blocks in flow charts should be understoodas representing specific aspects of module operation, or steps in theprocess. Alternate implementations are included within the scope of theexemplary embodiment of the present disclosure in which functions may beexecuted out of order from that shown or discussed, includingsubstantially concurrently, depending upon the functionality involved.

The broad teachings of the disclosure can be implemented in a variety offorms. Therefore, while this disclosure includes particular examples,the true scope of the disclosure should not be so limited since othermodifications will become apparent upon a study of the drawings, thespecification, and the following claims.

What is claimed is:
 1. A system comprising: a weighting function moduleconfigured to receive a plurality of samples of data, the data havingbeen obtained based on a head having read the data from a storage mediumwhile the head was at a given height above the storage medium, and applya weighting function to the plurality of samples to generate weightedsamples; a frequency sensor module configured to estimate, based on theweighted samples, (i) a first magnitude of a first frequency component,and (ii) a second magnitude of a second frequency component; and acontrol module configured to estimate, based on the first magnitude andthe second magnitude, the height of the head above the storage medium.2. The system of claim 1, wherein the control module is configured toadjust the height of the head above the storage medium in response tothe height of the head above the storage medium being estimated to be(i) greater than or equal to a first threshold or (ii) less than orequal to a second threshold.
 3. The system of claim 1, wherein theweighting function module is configured to generate the weighted samplesby multiplying the plurality of samples respectively by a plurality ofcoefficients of the weighting function.
 4. The system of claim 1,wherein the weighting function includes (i) a triangle window function,(ii) a flat-top window function, or (iii) a rectangular window function.5. The system of claim 1, wherein: the weighting function module isconfigured to apply the weighting function to the plurality of samplesin a time domain; and the frequency sensor module is configured toestimate (i) the first magnitude of the first frequency component, and(ii) the second magnitude of the second frequency component in afrequency domain.
 6. The system of claim 1, wherein the control moduleis configured to estimate the height of the head above the storagemedium based on a ratio of the first magnitude to the second magnitude.7. The system of claim 1, wherein: the first frequency componentcorresponds to a first harmonic, and the second frequency componentcorresponds to a third harmonic.
 8. The system of claim 1, wherein thecontrol module is configured to select frequencies of (i) the firstfrequency component, and (ii) the second frequency component.
 9. Thesystem of claim 1, wherein the frequency sensor module is configured toperform: a first discrete Fourier transform on the weighted samples toestimate the first magnitude of the first frequency component; and asecond discrete Fourier transform on the weighted samples to estimatethe second magnitude of the second frequency component.
 10. The systemof claim 1, wherein the data includes user data and not test data.
 11. Amethod comprising: receiving a plurality of samples of data, the datahaving been obtained based on a head having read the data from a storagemedium while the head was at a given height above the storage medium;applying a weighting function to the plurality of samples to generateweighted samples; estimating, based on the weighted samples, (i) a firstmagnitude of a first frequency component, and (ii) a second magnitude ofa second frequency component; and estimating, based on the firstmagnitude and the second magnitude, the height of the head above thestorage medium.
 12. The method of claim 11, further comprising adjustingthe height of the head above the storage medium in response to theheight of the head above the storage medium being estimated to be (i)greater than or equal to a first threshold or (ii) less than or equal toa second threshold.
 13. The method of claim 11, further comprisinggenerating the weighted samples by multiplying the plurality of samplesrespectively by a plurality of coefficients of the weighting function.14. The method of claim 11, wherein the weighting function includes (i)a triangle window function, (ii) a flat-top window function, or (iii) arectangular window function.
 15. The method of claim 11, wherein:applying the weighting function to the plurality of samples comprisesapplying the weighting function to the plurality of samples in a timedomain; and estimating (i) the first magnitude of the first frequencycomponent and (ii) the second magnitude of the second frequencycomponent comprises estimating, in a frequency domain, (i) the firstmagnitude of the first frequency component and (ii) the second magnitudeof the second frequency component.
 16. The method of claim 11, whereinestimating the height of the head above the storage medium comprisesestimating the height of the head above the storage medium based on aratio of the first magnitude to the second magnitude.
 17. The method ofclaim 11, wherein: the first frequency component corresponds to a firstharmonic, and the second frequency component corresponds to a thirdharmonic.
 18. The method of claim 11, further comprising selectingfrequencies of (i) the first frequency component, and (ii) the secondfrequency component.
 19. The method of claim 11, further comprising:performing a first discrete Fourier transform on the weighted samples toestimate the first magnitude of the first frequency component; andperforming a second discrete Fourier transform on the weighted samplesto estimate the second magnitude of the second frequency component. 20.The method of claim 11, wherein the data includes user data and not testdata.